function cr_output=cr_filter(input_sig, k)
    total=length(input_sig);
    cr_output(1)=0;
    for i=1:total-1
        cr_output(i+1)=(k*cr_output(i)+input_sig(i+1)-input_sig(i))/(k+1);    
    end
end